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1. Introduction 

One of the most interesting potential sources of low frequency gravitational waves 
(GWs) for LISA are the inspirals of stellar mass compact objects (white dwarfs, 
neutron stars or black holes) into supermassive black holes (SMBHs) in the centers 
of galaxies. The mass ratio is typically 10 : 10 6 , so these events are termed extreme 
mass ratio inspirals (EMRIs). Detection and parameter estimation for these events 
is likely to involve matched filtering. LISA will observe EMRIs for the last several 
years of inspiral prior to plunge, so the search templates will need to match the phase 
of the signal over a fewxlO 5 cycles. The extreme mass ratio ensures that templates 
of sufficient accuracy can be computed using black hole perturbation theory - 
the "self- force" formalism. 1 Such templates are not yet available, however, and will 
be very computationally expensive when they are ready. The number of templates 
required to cover the parameter space of possible EMRI signals is very large, 2 so 
there is a need for approximate models that are quick to generate while also being 
able to estimate the parameters of the source with sufficient accuracy that follow-up 
with more accurate waveforms is possible. 

One family of models are "adiabatic" templates, which are based on accurate 
evaluation of the dissipative part of the self-force, combined with the assumption 
that the orbital inspiral occurs slowly compared to the orbital period. 3,4 Adiabatic 
waveforms are likely to play a role at some stage of the LISA data analysis pipeline, 
but they are still computationally expensive. For scoping out LISA data analysis, 
waveforms must be generated in large numbers, so two families of approximate, 
quick-to-compute, "kludge" waveforms have also been developed. The "analytic 
kludge" is a phenomenological model based on Keplerian waveforms with relativistic 
inspiral and precession imposed. 5 The "numerical kludge" (NK) will be described 
here. These NK waveforms are sufficiently faithful that they may play a role in 
source detection for LISA and perhaps in source characterization. 



2. Numerical Kludge Waveform Model 

The NK family of waveforms are designed to be faithful models of true EMRI GW 
signals. The waveform parameters are the same as for true EMRI signals — using 
the NK model does not reduce the size of the parameter space or the number of 
waveforms required to cover it. For a given set of parameters, however, the model 
is much simpler to evaluate than a pcrturbative waveform and that is where the 
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computational savings for data analysis arise. Construction of an NK waveform 
is done in two stages — (i) generation of the trajectory for an object inspiralling 
through a sequence of quasi-geodesic orbits; (ii) construction of an approximate 
waveform for an object moving on this trajectory. 

2.1. Inspiral Trajectory Generation 

To construct the inspiral trajectory, we first compute the phase space evolution 
of the object, i.e., how the energy (E), angular momentum (L z ) and Carter con- 
stant (Q) of the orbit evolve with time. This is done by deriving suitable expres- 
sions for dE/dt etc. as a function of the orbital parameters, and then integrating 
them through phase space. The expressions we use are built on second order post- 
Newtonian (PN) results. 6 ' 7 Using 2PN results directly leads to pathological behavior 
for nearly circular orbits, but this can be corrected by amending the circular pieces 
of the fluxes. 7 The trajectories can be further improved by using fits to data derived 
from perturbation theory, i.e., based on solution of the Teukolsky equation. We have 
done this for circular orbits of arbitrary inclination, but not yet for generic orbits 
since perturbative data for such situations is only now becoming available. 4 The re- 
sulting phase space evolution equations are given in detail in Gair & Glampedakis 
2006. 7 For circular inclined inspirals, these fluxes agree with perturbative results to 
an accuracy of 1% for orbits with periapse greater than ~ 5M, and to an accuracy 
of < 5% for orbits with periapse greater than ~ 2M. For eccentric orbits, the fluxes 
agree to ~ 5% for orbits with periapse greater than 5M, but this increases to a few 
tens of percent for orbits that come very close to the central black hole. 

Once the phase space trajectory (E(t), L z (t), Q(t)) has been obtained, the in- 
spiral trajectory is derived by integrating the Kerr geodesic equations dr/dt = 
R(r,6,E,L z ,Q), dO/dt = 6(r, 9, E, L z , Q), dfi/dt = $(r, 6, E, L Z ,Q), with the 
time-dependent E, L z and Q inserted on the right hand side. We thus obtain the 
particle trajectory in Boyer-Lindquist coordinates, (r(t) , 0(t) , <f>(t)) . 

2.2. Waveform Construction 

After computing the particle trajectory in Boyer-Lindquist coordinates, we may 
construct a corresponding trajectory in a pseudo-flat space by identifying these co- 
ordinates with spherical polar coordinates. A waveform can then be obtained by 
supposing that there was a particle moving on such a trajectory in flat space, and 
using a weak-field GW emission formula. This approach is inconsistent in the sense 
that it neglects the stress-energy that is causing the particle to move on the tra- 
jectory, but it appears post facto to work well. We have constructed waveforms 
using the Press formula 8 (valid for weak-field, fast motion sources) and also us- 
ing the quadrupole and quadrupole-octupole formulae obtained by expanding the 
Press expression in v/c. Based on a balance between ease of computation and ac- 
curacy, it appears that the quadrupole-octupole formula is optimal. This waveform 
construction is described in more detail in Babak et al. 2006. 9 
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3. Application to LISA 

For both generic geodesic orbits, and for circular inclined inspiral orbits, the overlap 
between the NK waveforms and more accurate adiabatic waveforms is very high. 
For orbits with periapse greater than ~ 5M the overlaps are typically greater than 
95%, but this degrades for orbits that come deep into the strong field near the black 
hole. 9 The waveforms are sufficiently cheap to be generated in the large numbers 
required for LISA data analysis, while their high faithfulness suggests that they 
will also be able to constrain the source parameters quite well. NK waveforms are 
already being used for scoping out LISA data analysis, 2 and their high accuracy 
indicates that they could play an important role in source detection for LISA, and 
quite possibly for parameter estimation as the first stage of a hierarchical search. 

The NK waveforms can be further improved in several ways — (i) inclusion of 
PN conservative self- force corrections, i.e., the piece of the self- force that does not 
dissipate. We have already demonstrated how this can be done to lowest order for 
the simple case of circular inspirals in the Schwarzschild spacetime. 9 Inclusion of this 
effect will provide information on the relative influence of conservative corrections 
on the phasing of generic EMRI waveforms, currently a matter of some debate. 

(ii) Addition of "tail terms", i.e., the effect of radiation back-scattering off the 
background geometry. This can be done by expanding the Tcukolsky function and 
should help to improve the accuracy of the NK waveforms for strong-field orbits. 

(iii) Improvement of the flux expressions, i.e., dE/dt etc., for eccentric orbits by 
using fits to perturbative data. This will ensure the NK waveforms can match true 
EMRI signals for longer segments of the inspiral. These three improvements will be 
implemented in the future to further develop this model as a tool for data analysis. 
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